Rf-induced transport of Cooper pairs in superconducting single electron transistors in 

a dissipative environment 
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We investigate low-temperature and low-voltage-bias charge transport in a superconducting Al 
single electron transistor in a dissipating environment, realized as on-chip high-ohmic Cr microstrips. 
In our samples with relatively large charging energy values Ec > Ej, where Ej is the energy of 
the Josephson coupling, two transport mechanisms were found to be dominating, both based on 
discrete tunneling of individual Cooper pairs: Depending on the gate voltage V^, either sequential 
tunneling of pairs via the transistor island (in the open state of the transistor around the points 
Qg = CgVg = e mod (2e), where Cg is the gate capacitance) or their cotunneling through the 
transistor (for Qg away of these points) was found to prevail in the net current. As the open states 
of our transistors had been found to be unstable with respect to quasiparticle poisoning, high- 
frequency gate cycling (at / ~ 1 MHz) was applied to study the sequential tunneling mechanism. 
A simple model based on the master equation was found to be in a good agreement with the 
experimental data. 



INTRODUCTION 

During the last few years, the ultrasmall Josephson 
junctions have aroused the interest of researchers as con- 
venient objects for experimental investigations of macro- 
scopic quantum effects. A variety of their unique features 
formed the basis for possible applications in quantum 
information processing (see, e.g., review [l| and refer- 
ences therein). Cooper pair electrometry '2|, and metrol- 
ogy (accurate current sources on the basis of Cooper pair 
pumps 13). 

The Cooper pair dynamics in such junctions appears as 
a result of the competition between Coulomb and Joseph- 
son effects (see, e.g., 0), and it is strongly affected by 
interaction with the electrodynamic environment (for the 
review, see |^). In particular, the Cooper-pair transport 
is very sensitive to dissipative properties of the environ- 
ment. Quantitatively, these properties are characterized 
by the real part of the frequency-dependent impedance 
of the environment TieZ^uj). 

In the low-impedance case, |^(ti^)| ^ Rq, where Rq = 
/i/4e^ « 6.45 kf2 is the resistance quantum, the tunneling 
of Cooper pairs is elastic and occurs at zero bias voltage 
as a coherent process. The resulting net current presents 
the supercurrent as in the "classical" Josephson junc- 
tion because the Josephson phase across the junction 
is well defined. In most of the experimental systems, 
small but finite dissipation of the measuring circuitry 
|Z(u;)| ^ Rq causes phase fluctuations which transform 
the zero voltage supercurrent into a supercurrent peak of 
finite width located at small voltages (see Ref. [3| and the 
references therein). In contrast to the former two cases, 
tunneling of pairs in the junction embedded in a rather 
dissipative environment, KeZ^tu) ^ Rq, is related to an 
extensive energy exchange with the environment. Tun- 
neling of pairs in this case occurs incoherently and can 
be considered in terms of individual events. The first ob- 



servation of dissipative tunneling of single Cooper pairs 
was reported by Kuzmin et al. |3] for a simple system of 
a small single junction connected in series to a compact 
high-ohmic resistor. 

Several different transport mechanisms were reported 
for a two-junction superconducting device (single elec- 
tron transistor, SET) supplied by the high-ohmic resistor 
(see Fig.^l. One of the mechanisms - sequential tunnel- 
ing of pairs (ST) through either junction - was studied 
theoretically by Wilhelm et al. and experimentally 
by Kycia et al. ^10| and by Lu et al. These experi- 
ments were performed in an Al SET positioned on top of 
a two-dimensional electron gas with locally tunable dissi- 
pation (ground plane conductance) . Another mechanism 
based on a higher-order process of simultaneous tunnel- 
ing of pairs in both junctions - i.e., the cotunneling of 
pairs (CT) - was recently studied in detail Moreover, 
resonant tunneling of single Cooper pairs (RT) ITsI was 
found to be responsible for the charge transport at some- 
what larger bias voltages ~ 2Ec/e < 2Aai Here, 

2 

= Ec = 2^ is the charging energy, Cs = 2Ct + Cg the 
total capacitance of the transistor island, Ct the tunnel 
capacitance (in a symmetric device), Cg the gate capac- 
itance, and Aai ~ 200 /leV the superconducting energy 
gap of Al. The RT process involves resonant mixing of 
charge states across one junction, completed by tunnel- 
ing of Cooper pairs in the second junction with energy 
dissipation in the resistor. 

In our recent paper [l^ we had shown that at low 
temperatures and small voltage bias, CT dominates in 
the net current for gate charges Qg = CgVg, \Qg\ < e/2, 
where ST is blocked due to the Coulomb barrier > 2Ec, 
provided the Josephson coupling energy in the junctions 
i?j is small, <C Ec. The ST process becomes energet- 
ically favourable in the points around Qg ~ e mod (2e), 
giving rise to the "open" states of the SET. In the ex- 
periment, these states are often unstable towards spon- 
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FIG. 1: Equivalent circuit diagram of a symmetric supercon- 
ducting SET and different mechanisms of Cooper pair trans- 
port: sequential tunneling (ST), cotunneling (CT), and reso- 
nant tunneling (RT). 



taneous tunneling of non-equilibrium quasiparticles (the 
so-called, quasiparticle poisoning effect, see, e.g., and 
[l^'). Study of the crossover between ST and CT is there- 
fore a challenging task. 

In this paper, besides the observation of the CT cur- 
rent, we succeeded in realizing of the regime, where the 
ST process dominates over CT. We focus on the range 
of small bias voltages eV <C 2Aai, where tunneling of 
Cooper pairs dominates over that of quasiparticles, and 
RT is energetically unfavourable |l3j . In order to circum- 
vent quasiparticle poisoning of the island in the open 
states, we applied a high-frequency signal of sufficient 
amplitude to the gate and measured the time-averaged 
I—V dependencies. The measured curves were compared 
with the results provided by a simple model based on the 
master equation approach. We considered the current 
evolution along the gate cycle. As expected, it was char- 
acterized by the major input from ST in the open states. 
Our calculation also showed that the contribution of ST 
current, averaged over the whole cycle, dominates over 
that of CT. 



THEORETICAL BACKGROUND 

At voltages eV 2Aai ~ 400 ^eV, which are well 
below the thresholds for the different processes incorpo- 
rating tunneling of quasiparticles (see, e.g., Josephson- 
quasiparticle cycles in lljl 1 171 and jl8l|1. the tunneling 
of quasiparticles is rare |l9l |. and the net current is 
related to the transport of Cooper pairs. 

The Cooper pair dynamics in a superconducting SET 
with E] <C Ec and appreciable dissipation is qualita- 
tively similar to that of electrons in the normal-state SET 
described by the orthodox theory In particular, a 

quantum state of the whole system can be described by 



the number n of extra Cooper pairs on the island of the 
SET. The Josephson coupling in the junctions can be re- 
garded as a small perturbation in the total Hamiltonian 
of the system, giving rise to random tunneling events of 
different orders 0]. 

Since tunneling of Cooper pairs does not produce ex- 
citations, it is dissipating properties of the environment 
that define the functional dependence of a particular tun- 
neling rate on the electrostatic energy difference E as- 
sociated with the tunneling event. The corresponding 
probability function P{E) characterizes the ability of the 
environment to exchange an energy E with the tunnel- 
ing pair. For a given temperature, the shape of the P- 
function depends on the environmental impedance Z{uj) 
as well as on the tunnel capacitances of the junctions. 
In some important cases, P{E) can be expressed analyt- 
ically (see, for example, below). 

The process of the lowest order of perturbation is 
known as direct tunneling of a pair in either junction. 
It changes the charge state of the SET islan d by ±2e. A 
rate of direct tunneling can be expressed as 



1 by 
0: 



nE) = l^E^jPiE). 



(1) 



The validity range of the perturbative approach Eq. 
strongly depends on the dissipative properties of the en- 
vironment and is generally restricted to SETs with rel- 
atively weak Josephson coupling and/or large damping, 
i.e., when: 



EjPiE) < 1. 



(2) 



An analysis performed in Rcf. shows that in the 
case of zero temperature and a low-ohmic environment, 
ReZ(w) ^ i?Q, the probability P{E) grows infinitely 
around zero energies, and the condition Eq. |(2Jl is always 
violated at small voltages. On the contrary, at finite 
experimental temperatures T ~ 100 mK and/or in the 
case of moderate dissipation, ReZ(tj) ^ i?Q, the func- 
tion P{E) ~ ^E~^ is distributed over the broad range 
of energies so that the conditional equation jSJl can be 
reformulated as A = E'j /E^ ^ 4. For the junctions in a 
high-ohmic environment, ReZ(w) > i?Q, the conditional 
equation ^ is even less restrictive at low energies, where 
P is a power function of E (see Eq. (Ill) in Ref. j^). 
As a result, in the systems with noticeable dissipation, 
the approximation of discrete tunneling events can even 
be applied to junctions with relatively strong Josephson 
coupling. 

For small energies, the external impedance can be 
considered to be frequency-independent Z{uj) = Z{Q) 
(purely ohmic) and one obtains an analytic expression 
[22i | for the P-function for tunneling in individual junc- 
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which is vaHd for low temperatures and energies E <^ 
TioJz- The value of the threshold frequency uz of the 
ohmic approximation is defined for each experimental 
layout by all environmental components (including re- 
sistor, tunnel and stray capacitances and inductances ac- 
cording to the chip layout). In our experiment, the typi- 
cal value was evaluated ujz ^ 4x 10^° s~^ [hujz ~ 25 /xeV, 
cf. experimental section of this paper). In Eq. ||2Jl, T is 
a gamma function and 7 « 0.577 Euler's constant. The 
dimensionless resistance z = where p = Z(0)/i?Q, ac- 
counts for the factor 1/4 (in comparison to the case of 
a single junction or cotunneling) which arises from the 
partial decoupling effect of the tunnel junction from the 
environment due to the other junction (cf. network con- 
siderations in Ref. 0). The low-temperature approxima- 
tion requires fceT <C AEc/ttz which in most of the cases 
is valid in the whole mK temperature range. 

For tunneling in the left (right) junction of the sym- 
metric SET, the energy difference is: 



E{V, Q)^eV-4EJl± 



(4) 



it depends on voltage bias V and effective charge Q = Qg 
+ Qti where Qt = 2ne -I- me is a charge of n extra 
Cooper pairs and m extra quasiparticles on the island. 
Similar to the case of a normal-state SET (see, e. g., 
Ref. |2m), the energy relations in Eq. |0J can be mapped 
in the V — Vg plane, giving a 2e-periodic diamond-shaped 
stability diagram for Cooper pairs. As shown in Fig.|2Ia), 
this diagram is featured by a blockade domains with the 
stable number of extra pairs on the island, as well as by 
the regions of open states, where at least two alternating 
charge states are available for realization of sequential 
transport of Cooper pairs through the SET. At small 
biases, the values of the gate charge for the open states 
of the SET, e.g., Qg « e mod (2e) for m = 0, correspond 
to a small vicinity of the resonance points Q = ±e. 

Similar to the normal-state SET, it is also possible to 
calculate the net ST current at low temperatures, using 
the master equation approach. In particular, in the im- 
portant case of a small voltage bias, the stationary pop- 
ulation of the island is dominated by two charge states, 
n = 0, 1, with the occupation probabilities (t„ = (t„(i;, Ug) 
which can be found with the help of the detailed balance 
principle: 
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FIG. 2: (a) Stability diagram for Cooper pairs in a single- 
charge transistor for m = 0. The closed circles show the 
points of Cooper pair resonances in individual junctions, 
E — 0. The cycling amplitude Va is chosen to sweep over 
two open domains, (b) ST and CT currents, normalized to the 
Ambegaokar-Baratoff value I^^, were calculated as a function 
of the gate charge for the parameters of sample 2 (see Table^l 
and fixed V = 25 /iV and T = 100 mK. (c) Schematic repre- 
sentation of an rf-induced periodic sequence of current pulses 
resulting from cycling over two open domains (as shown in 
(a)). The amplitude of pulses given by the number of tunnel- 
ing pairs is fluctuating. 



and rj'^2(^) 



where Fs - F-(l)+r+(0)+F+(l)+F-(0) 
T^2{''^i ''^gi Qt) denote the rates in the left (right) junction 
in the positive (negative) direction and can be calculated 
with the help of Eqs. Q), (PJ and Q). The net current is 
obtained as a sum over the states n: 



I = -2e[aor+(0) -airr(l)] 
= -2e [a,r+{l)-aoT^{0)]. 



(6) 



An example of the ST current dependence on the 
gate voltages I{v,Vg), calculated for fixed bias voltage 
and the parameters of one of our samples, is shown in 
Fig. El;b). The net current, / - 0.14^^^ ^ 1 nA, where 
is the Ambegaokar-Baratoff critical current, is sig- 
nificant in the open domain between the resonant bound- 
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aries and decays exponentially with the height of the 
barrier \E\, see Eq. inside the blockade domains: 
/ cx exp(- /2A;bT). 

Simple analysis of the time-dependent master equation 
shows that the time r the system takes to reach its sta- 
tionary state, when moving between the operating points 
in the open domain (where the current is considerable), 
T ~ e/2I ^ 0.1 ns, is well below the experimental 

range of time ~ 10 ns of passing one open domain when 
the gate is cycled at / ~ 1 MHz. Making use of this 
fact, the quasi-stationary master-equation approach can 
be applied and the stationary ST currents considered for 
each point of the cycling trajectory, even when modeling 
high-frequency behaviour of the system. The gate cycling 
in particular will result in a regular pattern of ST current 
pulses, for example, 4 pulses per cycle which involves 4 
open domains (see Fig. |2Ic)), and the mean ST current, 
1st, can be obtained by averaging the values of I(v,Vg) 
over the cycling period. 

A contribution Id to the total net current due to 
higher-order processes, i.e., cotunneling of Cooper pairs, 
can be analyzed using the lines developed in Ref. It 
was shown in particular that an SET can be modeled by a 
single junction with an effective, gate- voltage-dependent, 
Josephson coupling Hamiltonian Htr{(t),Q), where is 
the overall Josephson phase. The coupling is 27r-periodic 
in 0, but its shape is generally different from —coscj). 
Expansion of the coupling Hamiltonian into Fourier se- 
ries (see Eq. (4) in Ref. 12] ) makes it possible to assess 
the individual contributions to the CT-currcnt from the 
different-order cotunneling channels, each one simulta- 
neously transporting a certain number k of Cooper pairs 
through the SET. Similar to direct tunneling, the rate 
formula in Eq. can be used to calculate the partial 
cotunneling rates due to each channel, with the substitu- 
tion E ~ 2keV , taking account of cotunneling of k pairs. 
The effective value of Josephson coupling energy for each 
channel is given by the corresponding Fourier coefRcient 
Ej^\q), and the dimensionless resistance k^p must be 
used in Eq. JSJ for the P-function. 

It was shown in Ref. that the lowest-order, i.e. 
single-pair, CT-process with k = 1 generally dominates 
in the CT-current through the SET, and the principal 
Fourier term can provide a reasonable approximation for 
the coupling Hamiltonian: Hti{<f>,Q) ~ —Ej{Q)cos{(f>). 
It can further be demonstrated that dropping-out of the 
higher terms causes an error of SI^ ~ 10 % at most which 
is obtained, in a symmetric SET, in the degeneracy points 
Q w ±e (see Eq. 6 in Ref. .12]). Around these points, 
the CT-current almost behaves like a Lorentzian peak, 
and at Q if decays more slowly than the exponen- 
tially suppressed ST current (see Fig. 12b)). As a result, 
the CT-current dominates in the net current inside the 
blockade domains. 



L= SEl EHT= 20.0 KU LJI}= L5 mm nflG= 'i 17.0 K PHDTD= 

2.DDum I 1 

sl_265_23_2/3inkin-ti lUO 




FIG. 3: SEM-micrograph of the Al/AlOx/Al-Cr - transistor. 
The configuration with two parallel Cr microstrips on each 
side was implemented for the purpose of directly measuring 
of Rcr as half the resistance between both right-hand arms of 
the device. 



EXPERIMENTAL SETUP 

The samples were fabricated by the shadow evapora- 
tion technique. The details of our fabrication process 
can be found in Each structure was built on three 
almost congruent metal layers (see Fig. ^ deposited at 
different angles through the same pattern in the mask 
with hanging bridges. The layers were of Cr, 7 nm thick 
and 0.1 /im wide as used for high-ohmic resistors, and 
two layers of Al, 30 and 35 nm in thickness, evaporated 
with intermediate oxidation of the first layer to form the 
tunnel barriers. The nominal tunnel junction area was 
40 nm X 80 nm. 

The results reported in this paper were obtained on 
two SETs (samples 2 and 3, see Table in a dissi- 
pating environment, realized as Cr microresistors with 
Z(Q) = i?cr ~ 3 and 12 kil as well as on a reference 
sample 1 without resistor, i.e., i?cr = 0. The correspond- 
ing values of the dimensionless resistance parameter were 
p « 0.5, 2 and 0. Due to the linearity and small lengths of 
the resistors it was possible to consider them as lumped 
elements (ohmic resistors) and neglect their distributed 
self-capacitance (« 65 aF//im, as earlier estimated for 
the similar layout ^2^). 

Ambegaokar-Baratoff critical currents of individual 
junctions were estimated, assuming symmetry of the 
transistors: I^^ = w 13—15 nA, where Rt is 

the tunnel resistance. The corresponding values of the 
Josephson energy were Ej = ~ 26 — 32 fieV, 

where — is the flux quantum. Roughly evaluated 
with the help of a similar device with i?cr = in the nor- 
mal state, the characteristic charging energy of the tran- 
sistor island was -Be ~ 150 to 170 ^eV. The ratio of the 
characteristic energies in all our samples, A « 0.15 — 0.2, 
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TABLE I: Sample parameters. 



Sample 


RT{kn) 


Rcrikn) 


(nA) 


E, (poV) 


-Ec (moV) 


-Ej/-Eo 


1 


20 





15 


32 


150 


0.22 


2 


20 


3 


15 


32 


150 


0.22 


3 


25 


12 


13 


26 


170 


0.15 



indicates that the SET charging energy dominated over 
the Josephson energy of the junctions, le-modulation 
period in the I —Vg curve without rf drive was w 20 mV 
yielding Cg « 8 aF. 

I — V characteristics were measured at the base tem- 
perature of the dilution refrigerator of about 10 mK. The 
bias and gate lines were equipped with microwave fil- 
ters made of Thermocoax® cables, 1.5m in length, each 
thermally anchored to the mixing chamber. The filters 
provided reliable electromagnetic isolation of the sample 
inside the shielded sample holder at frequencies above 
/ - 1 GHz 

EXPERIMENTAL RESULTS 

As discussed above, observation of ST current should 
be possible in the vicinity of the resonance points Q = 
±e. The parity effect in a superconducting SET (see, 
e. g., Refs. [13, and HI), basically makes it possible 
to reach these points through appropriately setting the 
gate voltage, provided quasiparticle poisoning (QPP) ef- 
fects are absent (the value of m is zero or, at least, fixed), 
and the experimental gate modulation is 2e-periodical. 
Unfortunately, the effect of QPP did not allow us to 
reach the resonance points and hindered the observation 
of ST currents in the dc regime: The observed periodicity 
was strictly le, corresponding to the range of the island 
charge — e/2 < Q < e/2, i. e., inside the ST blockade 
domain at low bias voltages. 

The effect of QPP is generally supposed to originate 
from various experimental non-idealities including sub- 
gap quasiparticle states [27| , out-of-equilibrium quasipar- 
ticles (28| , back influence of the SET electrometer on the 
island (box) j^^, etc. It involves uncontrolled excitation 
of quasiparticles in the outer electrodes and their tunnel- 
ing into the island of SET, thus changing m and mini- 
mizing the charging energy of the system. In particular, 
in the resonant points Qg — e mod (2e) and V — 0, the 
charging energy is minimized, when one extra quasipar- 
ticle enters the island resulting in Q = and Coulomb 
blockade of ST. The poisoning rate 7b depends on the 
number of excitations and the subgap tunnel resistance. 
Typically 7b - 16—1000 s'^ (see, e.g., Ref. 0). 

In Fig. ^ the dc I — V curves are shown for integer 
and half-integer values of Q in the system with consider- 
able dissipation (i?cr — 3 kf2, sample 2) in comparison to 
those for a similar SET fabricated without engineered dis- 
sipation {Rci = 0, sample 1). As expected for the block- 
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FIG. 4: Dc / — y curves of two SETs with and without on- 
chip resistor, samples 1 and 2, demonstrate the effects of both 
gate voltage and dissipative suppression of cotunneling. 



ade state, the dc currents were small (up to ~ 10 pA) and 
could be attributed to the higher-order CT process. In 
good agreement with the theory [l^l , the CT-current was 
modulated by the gate with a maximum-to-minimum ra- 
tio of «2 within the accessible range of— e/2 < Q < e/2. 

To avoid the effect of QPP and reach the open domain 
for ST at least for a short time, we applied a periodic gate 
signal of frequency / ~ 1 MHz which was well above a 
typical quasiparticle tunneling rate. Due to these fast 
oscillations, the actual values of Q along the gate sweep 
were expected to be less affected by the relatively rare 
quasiparticle tunneling than in the dc mode. As a result, 
the net current averaged over many cycles should include 
a substantial contribution due to ST in the open domains. 

The I ~ V curves measured in this regime are shown 
in Fig. |5l for sample 2 and several rf amplitudes of the 
gate signal Va, 1 < CgV^/e < 2. The most striking 
feature of these curves is the dramatic increase in the 
SET current, demonstrating a large contribution from 
the open domains involved in gate cycling (cf. scale of 
rf-induced current in Fig. [S] with that of dc / — 1^ curves 
of this sample in Fig. 0)). At fixed bias, the rf-induced 
current was a le-periodic function of the dc gate shift Vgo 
(cf. top inset in Fig. which shows that the effect of 
QPP was still substantial in a long time scale. 

A rough evaluation of the contribution of both ST and 
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FIG. 5: Rf-induced currents in sample 2 {Rcr ~ 3 ki}) for 
different cycling amplitudes. The cycle was centered in the 
middle of the blockade domain. Top inset: Rf-induced cur- 
rent vs. dc gate shift Vgo measured at 1/ = 50 fj,Y for cycling 
at / = 1 MHz and a relatively small amplitude Va ~ 20 mV, 
i. e., C^Va/g — 1. Bottom inset: I — V curve (symbols) at Va 
— 30 mV (CgVA/e = 1.5) compared with its fit Itot{V) (solid 
line) made at fixed Ec = 150 /ieV. The fitting values were: 
-Ej = 40 AtV, Rc, = 3.3 kQ., CgVA/e = 1.3, T =100 mK. The 
dashed line shows the calculated contribution of CT, IctiV). 



FIG. 6: Dependence of rf-induced current peak heights on the 
amplitude of rf cycling in sample 3 with significant dissipation 
[Rcr ~ 12 kO.). The gate cycling frequency was / = 3 MHz. 
The cycle was centered in the middle of the blockade domain. 
The shapes of rf-induced current peaks for this device are 
presented in the inset for different amplitudes Va = 12, 16, 
20, 40, 60, 72 mV (in the order of peak increase). Note the 
significant blockade at the smallest cycling amplitude Va = 
12 caused by insufficient gate sweep and small value of CT- 
current /ct ~ 1 pA at 1/ -100— 200/ieV. 



CT to the average current in the rf mode can be made 
with the help of duty-cycle weighting of typical open do- 
main values. For sample 2, these values were estimated 
to be 1st ~ 1-5 nA and let ~ 0.3 nA, respectively. Fur- 
ther assuming V = 2h and Va = 30 mV, the duty 
cycle is « 0.05 producing the estimates: 1st ~ 100 pA 
and let ~ 20 pA, respectively. The total current value 
^tot = 1st + let ~ 120 pA thus obtained agree with the 
order of magnitude of the experimental data in Fig. [S] 

For more accurate evaluation of 1st, we applied the 
quasi-static master equation approach (Eqs. lO)-®) and, 
for the total current Jtot, also took the cycle average let 
of CT current into account. We found satisfactory agree- 
ment between calculated and experimental I—V curves in 
the low voltage region, V < 25 fjY (see the bottom inset 
in Fig. [SJl. Noticeable discrepancy of the data at higher 
voltages in our opinion points to the upper limit of the 
frequency-independent behaviour of the environmental 
impedance assumed in Eq. Q for the calculation of the 
P-function. It is still difficult to explain the relatively low 
value of the cutoff frequency uiz — eV/fi ~ 4 x 10^" s^^ 
for our short (^0.5 /im - long) microstrips (cf. cut-off 
frequencies ~ 10^^ s^^ for ^eZ{uj) of similar microstrips 
in Ref. H). 

The height of the current peaks was found to quasi- 
periodically depend on the drive amplitude, with a ten- 
dency to saturate at higher amplitudes (cf . Fig. El show- 



ing the dependence for sample 3 with i?cr ~ 12 kVl). In 
our opinion, the observed quasi-periodicity on Va is a 
result of the step-wise involvement of multiple open do- 
mains in gate cycling, each contributing one pulse of ST 
current. In good agreement with this model, every pe- 
riod consists of two distinct segments: sharp increase in 
the peak height just above the values Va = (2n-l- l)e/Cg, 
n = 0,1,2,..., where an additional open domain becomes 
involved, and gradual reduction of the height along the 
other part of the period, corresponding to a shortening of 
the duty cycle with increasing amplitudes. The inset in 
Fig. El shows the rf-induced current peaks as a function 
of transport voltage V. The position of the peaks was 
evidently shifted to higher voltages as compared to the 
sample 2, which results from the shape of the P-function 
for a relatively large dissipation in sample 3. 



CONCLUSIONS 

We reported the experimental results on dissipative se- 
quential tunneling of Cooper pairs in SETs with compact 
on-chip resistors. Whereas Cooper pair cotunneling was 
the dominant dc transport mechanism for the gate volt- 
age range — e/2 < CgVg < e/2, we were able to observe 
rf-induced sequential tunneling of pairs realized with the 
help of high-frequency gate cycling resulted in a pulsed 
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ST current through the SET. We demonstrated that - 
in accordance with theoretical predictions - the integral 
contribution to the net current of ST exceeds that of the 
CT in the rf-cycling mode, indicating its dominance in 
the vicinity of the resonant points Q = ±e. By a simple 
model based on a quasi-static master equation approach 
it was possible to well describe the observed I — V curve 
at low voltages. 

Moreover, we demonstrated that, due to the effectively 
smaller environmental resistance seen by each single junc- 
tion of an SET ('^ g:-Rcr), the appreciable rf-induced ST 
currents can be realized even in systems with larger resis- 
tors of several i?Q , where the CT current is substantially 
suppressed. This effect, for example, is important for 
possible application of discrete Cooper pair tunneling to 
generate accurate currents by Cooper-pair pumps 

Another important result of our studies is the demon- 
stration of quasiparticle-poisoning-free performance of a 
superconducting SET at frequencies of / 1 MHz. This 
is important for the potential application of small Joseph- 
son junctions in quantum information processing, giving 
an estimate for a period of time of at least 1 /xs |30j 
as is available for reading out the final state of a charge 
qubit of the loop configuration with nonzero circu- 
lating supercurrent. 
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